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i The magnetic order resulting from a confinement-adapted Ruderman-Kittel-Kasuya-Yosida in- 

direct exchange between magnetic moments in the metallic phase of a AlAs/Gai-^Mn^As quantum 
well is studied by Monte Carlo simulation. This coupling mechanism involves magnetic moments 
and carriers (holes), both coming from the same Mn 2+ ions. It leads to a paramagnetic, a ferromag- 
netic, or a canted spin phase, depending on the carrier concentration, and on the magnetic layer 
width. It is shown that high transition temperatures may be obtained. 
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I. INTRODUCTION 
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During the last decade, due to the advances in the control of materials growth, and also in the techniques of 
characterization, a new interest arose in the study of the magnetic order in layered materials. This area is not 
restricted to magnetism in metals, but it also includes the study of magnetic semiconductor pseudo-binary alloys li 
Ai^ x M x B, where M stands, for a magnetic ion. These alloys are called Diluted Magnetic Semiconductors (DMS).Er 

Recently some groupsErEI succeeded in producing homogeneous samples of Gai-^Mn^As alloys with x up to 7% 
using low temperature (200 — 300° C) MBE techniques. Mn is a transition metal having its 3d level partially filled with 
■ five electrons, in such a way that it carries a magnetic moment of bfi/2, according to the Hund's rule. In the insulating 
phase, as in II- VI DMS, two Mn 2+ ions occupying the nearest neighbor positions are assumed to interact with each 
O ■ other via a super-exchange mechanism, resulting in an anti-ferromagnetic ordering of their magnetic moments. In the 
fee alloys, these interactions are frustrated, establishing the possibility of settling a spin-glass phase at low temperature. 
A double exchange mechanism which might stabilize a ferromagnetic coupling between the Mn ions in III-V DMS has 
been suggested by Akai,E3 but has not been confirmed in the EPR experiments performed by Szczytko et aljiil who 
did not observe trace of neutral Mn, concluding that the double exchange mechanism is not effective. 

The possibility of having a DMS based on GaAs opens a wide range of potential applications such as integrated 
magnetooptoelectronic devices. Besides its practical importance, this kind of DMS introduces an interesting problem: 
an Mn impurity in GaAs is an acceptor (it binds one hole), and at the same time it carries a localized magnetic moment. 
In the Gai-^Mn^As alloy Mn is, in fact, a strong p dopant,Bi3 the free hole concentration reaching even 10 20_21 cm -3 . 
At small Mn concentrations, the alloy is a paramagnetic insulator. As x increases it becomes ferromagnetic, going 
through a non-metal-to-metal transition for higher concentrations, and keeping its ferromagnetic phase. For x above 
7%, the alloy becomes a ferromagnetic insulator .0 In the metallic phase, depending on the value of x, the temperature 
of the ferromagnetic transition is observed in the range of 30-100 K, the highest values observed in DMS. The 
ferromagnetic order in the metallic phase is understood, at present, as resulting from the indirect exchange between 
the Mn 2+ ions due to the spin polarization of the hole gas. 

The aim of this work is to study the magnetic order resulting of the indirect exchange between magnetic moments. 
in a GaAlAs/Gai_ x Mn^As quantum well. A confinement-adapted Ruderman-Kittel-Kasuya-Yosida (RKKY)E^nL3 
mechanism is believed to be the most important interaction in such systems, if sufficiently strong doping is provided, 
as it is the case in metallic samples. It leads to an indirect exchange coupling between Mn 2+ ions, mediated by carriers 
(holes), which come from the same Mn 2+ ions. 

This article is organized as follows. In Sec. ||we present the calculation of the RKKY exchange for a confined Fermi 
gas in a semiconductor heterostructure. For the sake of relating our results with other previous ones, we explicitly 
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separate our calculations as intra-subband and inter-subband contributions. We emphasize that, in the quantum 
limit, i.e., when only the first subband is occupied, the intra-subband exchange is factoriaed into a purely 2-D RKKY 
exchange times a form factor determined by the architecture of the confining structures 



In Sec. Ill a Monte Carlo simulation is performed to determine the resulting magnetic phases and the relevant 
properties. Our calculations reveal that a ferromagnetic order may occur in a single DMS quantum well only beyond 
a minimum |W|idth of the magnetic layer, otherwise the sample is paramagnetic. This is in keeping with recent 
experimentsJ13 and is a consequence of the need of a certain number of magnetic neighbors before a ferromagnetic 
phase settles in. Depending on the well width and on the effective two-dimensional carrier concentration, a canted 
phase can occur, with a sizeable net low-temperature magnetization, < S > / S max , and a well behaved Edwards- 
Anderson order parameter q. The origin of the canted spin phase is investigated by analyzing the parallel and the 
perpendicular magnetizations, and the spin-spin correlation function components. The magnetic susceptibilities are 
calculated in the existing phases. 

Finally, in Sec. [V we summarize the results obtained, and comment about the expected magnetic order in the 
structures analyzed. 



II. RKKY INTERACTION IN A QUANTUM WELL 



The RKKY interaction between localized magnetic moments imbedded in a Fermi gas is a well understood problem, 
since its early developments almost fifty years ago. However, the new areas for experimental research brought into 
evidence some theoretical problems concerning the RKKY interaction in low dimensional systems, so far unexplored, 
such as purely 2-D and 1-D arrangements oi-i Magnetic moments, interaction between magnetic layers, magnetic 
moments in inhomogeneous electron gas, etcJlaEil The indirect exchange between localized magnetic moments in a 
quantum well mediated by a Fermi gas has been addressed several times. Basically, it deals with a confined electron 
(or hole) gas, therefore a quasi-two-dimensional system, being locally polarized by magnetic moments distributed in 
a layer. To our knowledge, Korenblit and Shender,EZl were the first to obtain a closed expression to the equivalent of 
the RKKY interaction in the limiting situation of a purely 2-D electron gas, although Kittellla obtained a numerical 
solution to this question earlier. Larsencj derived an expression for a general dimensionality, reproducing Korenblit. 
and Shender results for d=2. A-detailed calculation to obtain a closed expression in 2-D was shown by Beal-Monod.c3 
Gummich and da Cunha LimaEj studied the indirect exchange between magnetic impurities in a doped GaAs/AlAs 
quantum well in the diluted regime, obtaining a ferromagnetic interaction. Finally, another expression for a generic 
dimensionality has been derived by AristovH The approximation pf rthe Fermi gas in a quantum well by a purely 2-D 
system is seldom a reasonable choice. Helman and BaltenspergerOn3 treated the question of the polarization of an 
inhomogeneous electron gas in several circumstances, emphasizing the roles oi-jthe confined and extended states. The 
specific case of a DMS quantum well was addressed recently by Dietl et al,E3 but they assumed that the magnetic 
moments are spread all over the region allowed to the carriers, and in that case the inter-subband contributions to 
the Curie- Weiss temperature cancel out in a mean field approximation. 

The interaction potential between a Fermi gas and a set of localized magnetic moments at positions Ri is well 
described by the Hund-type exchange potential: 

ffcx = s(r)S(r- Ri), (1) 

i 

where Si is the spin of the magnetic moment at position Ri, which will be treated as a classical variable, and s(r) 
is the spin operator of the fcrmion at r. I is the sp — d interaction.^ If ipa(r) and "4>\{f) describe the fermion field 
operator for spin cr, then 

<f{?) = \(^(r)^(r) ~ ^(r)^i(f)), (2) 



s +(f)=^(r)^(r), (3) 

S -(f)=ViJ(r)^ T (f) ! (4) 

with the usual definitions of s + = s x + is y , and s~ = s x — is y . Instead of free fermions in a 3-D space, the electrons 
and holes in a ssemiconductor heterostructure are confined in the growth direction, assumed to be the z-axis, due to 
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the mismatch of the conduction and valence band edges. Since they are free particles in the plane perpendicular to 
that growth direction, i.e., in the plane parallel to the semiconductor interfaces, their field operator are given by: 



n.k 



(5) 



where A is the normalization area, A: is a wavevector in the plane (x,y), rj is the spin tensor, <f> n (z) is the envelope 
function which describes the motion of the fermion in the z-direction, and c n g is the fermion annihilation operator 

for the state (n,k,a). R represents a vector in the 2-D coordinates plane (x,y). The usual RKKY perturbation 
calculation up to second order leads to the correction on the ground state energy of the system formed by the set of 
(classical) localized moments and the Fermi gas : 



where 



6E& = 6E™ + 5E™ , 



(6) 
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Eqs. (0) and (||) are, respectively, the self-energy term and the RKKY exchange in the form they assume for confined 



fcrmions. The coordinates (Ri,Zi) describe the position of the impurity i in the plane (x,y). 
direction; q is a two-dimensional wavevector. x n ' n io) is the equivalent to the Lindhard function 
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Eqs. (0) and (||) are used to define the exchange Hamiltonian: 

H ex — ^ ^ JijSi ■ Sj 



(10) 
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A. Intra-subband terms 

To our knowledge, complete calculations of Eq. (|^) have only been performed for intra-subband transitions, us- 
ing different approaches. For the sake of completeness, we will show how this is achieved in our treatment. The 
contribution of a subband n to the exchange reads: 

^ = (2l) 2 1 MZl) |2 ' MZl) I' E 2cos (« ' 4)x"' n (9l- (12) 

q 

We observe that, in the so-called quantum limit, when only the first subband (n = 0) is occupied, the difference 
between Eq. (|l^) and the indirect exchange mediated by a 2D electron gas comes from the non-uniform charge 
density in ther-aonfining direction z. Actually, in that case, Eq. ( |l2| ) factorizes into a form factor T 1 ^ and a purely 
2-D exchange jla 
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Jg ) =^jg- D) . (13) 

where 

=| M«) l 2 l M*j) I 2 (14) 

and 

It is easy to show, by using the dimensionless variables x — kRij , and y = gRg, that the Fourier transform of the 
modified Lindhard function, appearing in the summation in q at the rhs of Eq. (|I2|), becomes: 

4mM 2 f°° r^-Rij r*/ 2 \ 

x n (Rv) = , t2n , / dyyMv) / dxx / — — ti> ( 16 ) 



4^ = E 2cos (^4)x n ' n ®. (15) 



TT 3 Ti 2 R?j Jo yyuyy 'J Jo V-4x 2 cos 2 

where x"(Rij) — J2q^ cos i ( l ' Rij)x n ' n {<t)- The transversal effective mass, m|, is assumed as isotropic in the plane 

parallel to the interfaces. As usual, kp — y/2m*(Ep — e n )/%. Performing the 4> integral and changing variables 
again (y/2x -> y), 

2m* A 2 r ki F ] - R ii r°° i 
X n (^i) = "Sbf / ^ / dy^o(2xy)-==. (17) 

7r 2 n i? 2 - Jo 7i Jy 2 - 1 



The integral on y is straightforward S 



X n (%) = -— 1— ■ / &xJ (z)Af (z). (18) 

7T?l if?- 



m\A 2 



After performing the integral on x, Eq. ( Jlq ) results in: 

x n (%) = - ! rr4" )2 ^(4 l, ^)A f o(4 n) ^) + Ji(4 n) ^j)^i(4 n) «ij)]- (1 9 ) 

This piBfipsian for the real space Lindhard function has been derived in a different context, by several 
authors .EllESEjly The final expression for the intra-subband exchange becomes: 

^ = -(0 2 ^4 n)2 i^) i 2 i^-)i 2 x 

[ Jo (fc£° JZy ) JVb (4 n) % ) + * (4 n) % )^i (*&° ^ )] • ( 2 °) 



B. Inter-subband terms 



The contribution of the inter-subband terms cannot be expressed easily in a closed form. Starting over from Eq. 
(|ll|), and using the same approach as in Ref. [13], we arrive to: 

' / / \ 2 1 f°° 

J « UJ -^M*)M*i)#(*i)M*i)]y dqqF n , n ,(q)Jo(qRij), (21) 

where we used 

f „, ( „ . A _£±A^_ fl(w _ Ef) , (22) 

and A„/ jn = 2toj • — E n )/fi 2 . The integral in Eq. (p2|) is, then, straightforward: 



Fn,n'(Q) = (1 - - V 1 " ( 2 2 5l )g )^^ 2 + A «'.« ~ W^Wn' ~ E F ). (23) 

27T?l <T V + A "',n 
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III. MONTE CARLO SIMULATION: MAGNETIC ORDERING 



In order to determine the possible magnetic order in GaAs:Mn quantum wells, we have performed extensive Monte 
Carlo simulations. Classical spins Si, randomly distributed on the cation sites with concentration x, are assumed to 
interact through the RKKY exchange Hamiltonian defined by Eq. (JTc|) . 

In the present work we have focused our attention on metallic single quantum wells with the Mn concentration 
x = 5%, and we have neglected possible (anti-)ferromagnetic interaction between the nearest neighbors and the next 
nearest neighbors pairs. The RKKY exchange interaction derived in Sec. || is assumed to be effective within a cutoff 
radius which we have taken as R c = 4a, R c = 2a, and R c = a, where a is the fee lattice parameter of GaAs. This 
makes the smallest value assumed for R c nearly equal to the hole mean free path estimated from bulk transport 
measurements!] The highest value, R c = 2.2nm, amounts to 3-4 values of that mean free path. The consequences of 
the cutoff radius on the results will be discussed below. 

The calculation is performed in a finite box, whose axes are parallel to the [100] directions. Its dimensions are 
L x = L yi and L z — Na/2, and N is the number of DMS monolayers (ML) in the barrier. Periodic boundary conditions 
are imposed in the (x, y) plane. The lateral dimensions are adjusted in such a way that the total number iV s of spins 
is about 4400, for all samples with different L z . The initial spin orientations are randomly assigned. At a given 
temperature, the energy of the system due to RKKY interaction is calculated, and the equilibrium state for. a given 
temperature is sought by changing the individual spin orientation according to the Metropolis algorithmic A slow 
cooling stepwise process is accomplished making sure that the thermal equilibrium is reached at every temperature. 
The resulting spin configuration is taken as the starting configuration for the next step at a lower temperature. 

For everjz-temperature, the average magnetization < M > and the Edwards- Anderson (EA) order parameter q are 
calculated. L3 The latter is defined as 



2 = 1 \ Q t'—tQ J 



1/2 

(24) 



where a — x, y and z. In order to avoid spurious results in obtaining the average over a large time interval t, a 
summation on t' is performed starting from a time to, when the system already reached the thermal equilibrium. 

In our calculations we used the value Nq/3 = — 1.2 eV, taken from Ref. |3Jj, for which we obtain transition tem- 
peratures injjpod agreement with the experimental data. Recently, the value No/3 = —0.9 eV has been obtained 
theoreticallyCJ confirming the result of Ref. |3^. The earlier estimate!] \Nq(3\ = 3.3 eV is probably too high. 

Monte Carlo calculations have been performed in sixteen samples, as shown in Table ||, the well widths varying 
from 25 A to 100 A (9 ML to 35 ML). In sample #01, for a well width of 50 A and assuming R c = 8 ML, we tested 
the effect of a small hole concentration, making p just 1% of x. This amounts to p » 1.1 x 10 19 cm~ 3 . We found that 
the spins can be arranged in a ferromagnetic phase even at this carrier concentration. For that sample the calculation 
gives a transition temperature near 27 K. 

In Fig. ^ the normalized average magnetization is shown for samples #02 to #07 with different well widths, but 
with a fixed carrier concentration of 10% of x, i.e., p w 1.1 x 10 20 cm~ 3 , and the same cutoff R c = 8 ML. For a very 
thin well (25 A), we found that the sample is paramagnetic. The EA order parameter for these samples is shown in 
Fig. ||. It can be observed that, for sample #02, there is no phase transition. All the other curves in Fig. || (samples 
#03 to #07) are characteristic of an ordered phase. Raising the well width (starting from 35 A), the number of 
interacting neighbors increases, and the sample shows successively a ferromagnetic phase (samples #03, #04, and 
#06) and a canted spin arrangement (samples #5 and #7). The value chosen for the cutoff radius is larger than the 
first zero of the J^ KKY , so antiferromagnetic interactions are turned on. With these choices of p and i? c , depending 
on the well width, the antiferromagnetic interactions can settle a fraction of the Mn magnetic moments antiparallel. 
This is the origin of the canted spins phase. The final average magnetizations in Fig. |l| at T = K are only a fraction 
of the maximum magnetization, around 60% for sample #5, and 70% for sample #7. The transition temperatures 
were found in the range of 40 to 50 K. No spin glass phase with vanishing magnetization was found. 

The existence of the canted phase requires a careful analysis of a possible dependence on the cooling process 
initial conditions. In order to clarify that point, we performed two additional simulations on sample #5, starting 
at T = 25 K, where the sample shows already a significant partial alignment of spins, and we proceeded with the 
slow cooling down process. In the first simulation, we assumed a starting configuration in which all spins are aligned 
perpendicularly to the interface. In the second, the spins alignment is made parallel to the interfaces. This choice 
of a rather low starting temperature for cooling is necessary, otherwise the thermal excitation would immediately 
randomize the initial configuration. The results are shown in the inset of Fig. [j]. We observe that the appearance 
of the canted phase does not depend on the choice of the starting configuration, and the three simulations converge, 
within statistical fluctuations, to the same value of the magnetization at every temperature step. 
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In Fig. |3j the magnetization as a function of temperature is shown for samples #08 to #12, with a higher carrier 
concentration, p — 0.25x, but keeping the same cutoff radius of 8 ML. The EA order parameter for these samples 
shown in Fig. | as a function of temperature, gives evidence for the existence of ordered phases. The Fermi wave 
number increases with carrier concentration, what decreases the in-plane distance corresponding to the first zero of 
the RKKY interaction. In consequence, in some samples, the magnetic moments order in the canted spin phase 
(samples #10 to #12). In other samples, however, the ferromagnetic interaction prevails and the total magnetization 
is reached (samples #08 and #09). Notice that the canted phase appears here already for L = 50 A, while the phase 
is still ferromagnetic for that width when p = O.lx. The transition temperatures were estimated to lie between 30 K 
and 50 K. 

In Fig. H we explored the effect of the cutoff radius on the spin ordering with p = 0.25a;. A cutoff radius R c — 2 ML 
was used in samples #13 (L = 60A) and #14 (L = 80A), while R c = 4 ML was used in samples #15 (L — 60A) and 
#16 (L — 80A). The former is too small, resulting in the fact that no net magnetization is allowed. The respective 
EA order parameters (Fig. ||) are typical of a paramagnetic phases for samples #13 and #14. The choice of a larger 
R c ordered spins in a ferromagnetic phase, with transition temperatures calculated to be 35 K for the sample #15, 
and 80 K for sample #16. The cutoff radius of 4 ML is smaller than the first zero of J^ KKY . In this situation the 
canted spin arrangement is not allowed, and the sample can be either paramagnetic or ferromagnetic. 

Finally in Fig. [7] the magnetic susceptibility for sample #07 is presented, calculated from the equilibrium magne- 
tization fluctuations. Notice that the peak in the in the susceptibility indicates the same T c as estimated from the 
magnetizations and EA order parameter curves (Figs. |l| and || respectively). 

IV. DISCUSSIONS AND FINAL COMMENTS 

The results of the Monte Carlo simulations indicate that, besides the choice of the range of the interaction (the 
cutoff radius R c ), two parameters are determining in the magnetic ordering in these heterostructures: the magnetic 
layer width L, and the carrier concentration p. r-.r—. 

It is observed from resistivity measurements, that the holes have a small mean free pathQH in these materials. The 
criteria for choosing the cutoff in a Monte Carlo simulation must take into account natural scales of the interaction. 
These scales are the transport mean free path (since the RKKY interaction is based on the very existence of free 
carriers), and also the spin coherence length. We tested different i? c 's, in the range of the mean free path estimated 
from resistivity data. The influence of this parameter on the magnetic order is simple. If R c is smaller than the first 
zero of J^j KKY (a proper choice for the case in which the transport mean free path or the spin coherence length are 
small), there are two possibilities of magnetic order: ferromagnetic or paramagnetic. For larger i? c , on the other hand, 
corresponding to the cases where both the transport mean free path and the coherence length are large, a canted 
magnetization may be observed. 

In all the explored samples, no spin glass phase was found. This is presumably due to the fact that while the spin 
frustration exists, as witnessed by the occurrence of canted spin phases, it is not strong enough to produce a spin 
glass phase, as in canonical metallic spin glasses. A spin glass phase in these DMS structures would probably require 
a much higher carrier concentration. 

In what concerns the influence of the width of the quantum well, we conclude that, for thin layers, the number 
of interacting ions is small within the cutoff radius, and the sample is paramagnetic. When L becomes larger, 
the number of interacting ions increases and a collective magnetic ordering may be observed. The fact that the 
appearance of a magnetic order occurs only above a minimum thickness of the magnetic layer has already been 
observed experimentallyO 

Since the RKKY interaction oscillates with the argument (kpR), which depends on the carrier concentration, raising 
p produces a change on fcp, increasing the number of oscillations of J^ KKY . Therefore, antiferromagnetic interactions 
can be turned on, resulting in all kind of couplings. In this situation, ferromagnetic and antiferromagnetic interactions 
compete in establishing the magnetic order, which depends on the other sample characteristics, resulting in a total 
or in a partial alignment of the Mn magnetic moments. The occurrence of partial magetization (about 40%) has also 
been observed in samples of (In,Mn)As/(Ga,Al)Sb.EJ 

To conclude, we believe that the RKKY mechanism explains the high transition temperatures experimentally 
observed in Gai-^Mn^As heterostructures, at least in the metallic phase. Additionally, it explains, as it becomes clear 
after these Monte Carlo simulations, the occurrence of samples showing a partial magnetization at low temperatures. 
The possibility of having a ferromagnetic phase in samples with a low Mn concentration, i.e., in the ferromagnetic 
insulator Gai-^Mn^As, remains to be explained. 
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TABLE I. Sample characteristics: L is the well width; r is the ratio of the carrier concentration to the Mn concentration; T c 
is the transition temperature for the phases: F: ferromagnetic, P: paramagnetic, C: canted spin; R c is the cutoff radius of the 
RKKY interaction, in number of monolayers. 



sample 


L(A) 


r 


Rc 


phase 


Tc (K) 


#01 


50 


0.01 


8 


F 


27 


#02 


25 


0.10 


8 


P 


< 1 


#03 


35 


0.10 


8 


F 


50 


#04 


50 


0.10 


8 


F 


50 


#05 


60 


0.10 


8 


C 


50 


#06 


80 


0.10 


8 


F 


50 


#07 


100 


0.10 


8 


c 


40 


#08 


25 


0.25 


8 


F 


50 


#09 


35 


0.25 


8 


F 


50 


#10 


50 


0.25 


8 


C 


50 


#11 


60 


0.25 


8 


c 


40 


#12 


100 


0.25 


8 


c 


30 


#13 


60 


0.25 


2 


p 


< 1 


#14 


80 


0.25 


2 


p 


< 1 


#15 


60 


0.25 


4 


F 


35 


#16 


80 


0.25 


4 


F 


80 



FIG. 1. Normalized magnetization vs temperature for samples #02 to #07 indicated in Table IE The inset shows the low 
temperature magnetization with three starting spin configurations: random, parallel and normal to the layers. 



FIG. 2. EA parameter vs temperature for samples #02 to #07 indicated in Table 



FIG. 3. Normalized magnetization vs temperature for samples #08 to #12 indicated in Table 



FIG. 4. EA parameter vs temperature for samples #08 to #12 indicated in Table 



FIG. 5. Normalized magnetization vs temperature for samples #13 to #16 indicated in Table 



FIG. 6. EA parameter vs temperature for samples #13 to #16 indicated in Table |. 

FIG. 7. Magnetic susceptibility vs temperature, calculated from equilibrium magnetization fluctuations, for sample #07 
indicated in Table M. 
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